Coulomb-blockade-induced bound quasiparticle states in a double-island qubit 
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We study the low temperature thermodynamic properties of a superconducting double-island 
qubit. For an odd number of electrons in the device, the ground state corresponds to the intrinsic 
quasiparticle bound to the tunneling contact. The ground state is separated from the continuum 
of excited states by a finite gap of order of the Josephson energy E.j. The presence of the bound 
quasiparticle state results in a nonmonotonic temperature dependence of the width of the tran- 
sition region between Coulomb blockade plateaus. The minimum width occurs at the ionization 
temperature of the bound state, Ti ~ Ej / ln(V AEj / 8) , with A and 8 being respectively the su- 
perconducting gap and single particle mean level spacing in the island. For an even number of 
electrons in the system, we show that the Coulomb enhancement of the Josephson energy can be 
significantly stronger than that in the case of a single grain coupled to a superconducting lead. If 
the electrostatic energy favors a single broken Cooper pair, the resulting quasiparticles are bound 
to the contact with an energy that is exponentially small in the inverse dimensionless conductance. 

PACS numbers: 73.23.Hk, 74.50.+r, 85.35.Gv 



Some of the most promising qubit realizations are 
based on superconducting nanodevices 0, 0, 010, IE • 
These include charge El- [J El Q and phase @ qubits. 
A charge qubit consists of a small superconducting island 
coupled by a tunneling contact either to another super- 
conducting grain or a bulk superconducting lead. At low 
temperatures quantum fluctuations of the island charge 
strongly influence the properties of such devices. 

We study the effect of quantum charge fluctuations 
on the thermodynamic properties of a double-island 
qubit 0, H. The device is schematically depicted in 
Fig. ^ It consists of two superconducting islands cou- 
pled to each other by a tunnelling contact and capac- 
itively coupled to external gates. In a typical experi- 
mental realization the probability of electron tunneling 
between the qubit and the surrounding environment on 
the experimental time scale is negligible. Accordingly we 
consider the total number of electrons in the device to 
be fixed Q • This number may either be even or odd and 
can be randomly changed by heating the device. 

In the absence of tunnelling between the grains the 




FIG. 1: The two islands of the qubit are coupled to the gates 
which are maintained at the voltages V g i and V g 2- 



electron number on each island is quantized and changes 
abruptly at the charge degeneracy lines in the V g \-V g i 
plane. Depending on the values of the pairing gap and 
the charging energy electrons can be transferred between 
in the islands either one by one or in pairs. At finite tun- 
nelling between the grains quantum charge fluctuations 
lead to rounding of Coulomb blockade steps in the volt- 
age dependence of the island ch arg e. This phenomenon 
was studied theoretically [1 [1 Ell El It2| and experi- 
mentally EH for different single electron devices. 

In this paper we determine the gate voltage depen- 
dence of the island charge for a double- island qubit. We 
also study the Coulomb enhancement of the energy of 
Josephson coupling between the grains. For a fixed total 
number of electrons in the qubit the Hamiltonian of the 
system can be written as 
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where t^ p are the tunneling matrix elements, a = ±1 is 
the spin index, and a p /k,<r is the annihilation operator for 
an electron in state k, a in the left island or state p, a in 
the right island respectively. The operator of the electron 
number difference (in units of electron charge) between 
the left and right islands is denoted by n, Ec is the corre- 
sponding charging energy, and the dimensionless param- 
eter N is given by the appropriate linear combination of 
the voltages on the left and right gates. Finally, H de- 
notes the Hamiltonian of the two grains in the absence 
of tunneling and charging interaction. For simplicity we 
consider Ho to be of the BCS form and assume that both 
islands have identical pairing gaps, A, and single particle 
mean level spacings, 8. 

We assume that (\tk P \ 2 )kp/S 2 <C 1 , where (...)k P 
denotes averaging over the Fermi surfaces of both is- 



lands. Then the dimensionless conductance 
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8tt 2 ^ tkp s [ " >kp of the tunneling contact, with R being its 
normal state resistance, obeys the condition g <C 8ir 2 . 
Further, we restrict our analysis to low temperatures, 
T < T* = A/ln(A/<5), where no thermal quasiparticles 
are present |l4|. 

First we consider the case of an odd number of elec- 
trons in the device. To be specific we assume it to be 
equal to unity and study the gate voltage dependence 
of the average island charge difference (h) for \N\ <C 1. 
In this regime the charge states with n = ±1 are nearly 
degenerate and quantum charge fluctuations are strong. 

In the absence of tunneling the eigenstates of the 
Hamiltonian, Eq. (yi are characterized by n and the mo- 
menta and spins of the quasiparticles present in the is- 
lands. For an odd number of electrons there is precisely 
one quasiparticle in the qubit at T <T*. We denote the 
states with n = 1 and the quasiparticle in state k, a in 
the left island by \k, a; 0) and those with n = — 1 and the 
quasiparticle in state p, a in the right island by \0;p, a). 

For weak tunneling the most straightforward way to 
proceed is to use perturbation theory. In second order 
in tunneling scattering matrix elements between states 
\0;p,a) and \0;p',a), and fc,cr;0) and |fc',<r;0), appear. 
Their magnitude ~ iihsl ^ fc p ^J~Z^kZ significantly exceeds 

the energy distance between such states, ~ 6 2 /A, even 
far from the charge degeneracy point. Thus in contrast 
to the case of a normal grain y| virtual tunneling pro- 
cesses strongly mix different quasiparticle states even for 
{\tkp\ 2 )k P /5 2 <C 1. For |iV| -c 1 only the two nearly de- 
generate charge states need be considered. In the spirit 
of degenerate perturbation theory we look for the 
ground state wave function in the form of a linear com- 
bination of low energy states 
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We rewrite the tunneling Hamiltonian (|lc|l in terms 
of the quasiparticle creation and annihilation operators 
with the aid of the Bogolubov transformation a p tk),a = 

Up(k) a p(k),o- + aV p(k) a -p(k)-a ( here U p(k) End V p{k) are 

the standard BCS coherence factors) and project it onto 
the subspace of low energy states given by Eq. © . The 
energy spectrum of the resulting problem can be deter- 
mined exactly. The density of states in the continuum is 
unaffected by the presence of tunneling. In addition to 
the continuum states one bound state below the contin- 
uum appears. Though its wave function, Eq. (J2J), depends 
on the specific form of the tunneling matrix elements tk P 
its energy in the weak tunneling approximation can be 
expressed in terms of the dimensionless conductance of 
the tunneling contact, g. We obtain for the ground state 
energy 
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where Ej — gA/8 is the Ambegaokar-Baratoff ex- 
pression for the Josephson energy. The binding energy, 
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FIG. 2: The dependence of the island charge difference on the 
rescaled gate voltage 2EcN / Ej for temperatures T = (a), 
T = OMEj (b), T = 0.2Ej (c), T = Ej (d). The tempera- 
ture dependence of the step width is clearly non-monotonic. 



i.e. the gap between the bound state energy and the 
bottom of the continuum is given by 
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where we introduced a rescaled gate voltage N = N. 

Using the obtained energy spectrum we calculate the 
free energy F and find the gate voltage dependence of the 
island charge difference n(N,T) = (h) = —je^TiTj + 
where (. . .) denotes the thermodynamic average. The 
results are conveniently expressed in terms of the rescaled 
temperature T = T/Ej, binding energy 1b = £b/Ej and 
gate voltage N. Denoting the zero temperature result by 
n(N,0)=n (N), 



n (N) = 
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we can express the island charge difference as 
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Here we introduced the notation a 

The gate voltage dependence of the island charge is 
plotted in Fig.|21for several temperatures and typical ex- 
perimental parameters 4j, A w 2K : g w 1, S f» 10~ 4 K. 
Remarkably, the temperature dependence of the step 
width is non-monotonic and shows a minimum at the ion- 
ization temperature of the bound state Ti ~ Ej /a. Qual- 
itatively, the non-monotonic temperature dependence of 
the step width can be understood as follows. At T < Tj 
the quasiparticle resides predominantly in the bound 
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state, and the width of the step is determined by the 
binding energy, SN ~ EjjE c . At T > 7} the quasiparti- 
cles resides in the continuum, and the width of the step 
is determined by the temperature, SN ~ T/Eq- Since 
according to the Saha rule the ionization temperature Tj 
is smaller than the binding energy Ej the temperature 
dependence of the step width is nonmonotonic. 

We now turn to the case of an even number of electrons 
in the device, which we take to be 2 to specific. For 
brevity we set T = from here on. 

First we study the dependence of the strength of the 
Josephson coupling between the grains on A — 2Ec for 
A > 2Ec- We consider the vicinity of the degeneracy 
point (\N\ <C 1) between the two charge states, n = ±2, 
with one extra Cooper pair in the system. For A > 2Ec 
the energies of these two states lie below the bottom of 
the continuum of states with n = and one quasiparticle 
in each island, see Fig [3] 

For A— 2Eq sufficiently large the Josephson energy can 
be obtained from perturbation theory in tunneling 
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where e Hp) = 

At A — > 2Ec the integral in Eq. Q diverges logarith- 
mically. This divergence arises from tunneling through 
the states of the continuum with £ fc ( p ) <C A. These states 
become strongly hybridized with the two discrete states 
with n = ±2. To find the Josephson energy with log- 
arithmic accuracy we look for the wave function in the 
form of a linear combination of the low energy states, 

|*) = * 20 |2,0) + * 02 |0,2)+^ C Plkl(T \k,a;p,-a}. (8) 

Here |2,0) and |0,2) denote the states with n = 2 and 
n = — 2 respectively, and \k, a;p, —a) denotes the n = 
state with quasiparticles in state k, a in the left dot and 
p, —a in the right one. Solving the Shroedinger equa- 
tion with this wavefunction is straightforward. The con- 
tunuum spectrum is unaffected by the presence of tun- 
neling. The energies E of discrete spectrum can be found 
from the equation 
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Here En = EcN 2 is the electrostatic energy of the 
fc, <r;p, —a) states, and i?20(02) = Ec{2 T -^) 2 are those 
of the |2,0) and [0, 2) states respectively. For A > 2Ec 
this equation has two solutions. The Josephson energy 
Ej is given by the energy splitting between these states 
at the degeneracy point N = 0, 



Ej = -Ej In 
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where Ej is the Ambegaokar-Baratoff expression for the 
Josephson energy defined below Eq. ©. 




FIG. 3: Energy diagram for the Josephson energy calculation. 
Electrostatic energies of the relevant charge states are plotted 
as a function of N. The states with n = ±2 are coupled by 
tunneling via the continuum of states with n = 0. 



For Ej < A — 2Eq <§C A this expression reproduces 
the logarithmic asymptotic of the perturbative result J7J. 
For < A — 2Eq < Ej perturbation theory breaks down 
and the Coulomb enhancement factor of the Josephson 
energy saturates for small g at the value 
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Equation (|10|l shows that the Coulomb enhancement 
factor of the Josephson energy in a double-island qubit 
can be much larger than that for a single grain connected 
to a bulk lead [ItI Il8|. which cannot exceed 1.3. This 
agrees with the recent experimental observation of an en- 
hancement by a factor of ~ 3 in double- island qubits Q . 

For A — 2Ec < the discrete level with the higher 
energy disappears as a result of merging with the contin- 
uum, whereas the lower level continues to be separated 
by a finite gap from the continuum and represents the 
ground state of the system. The gate voltage dependence 
of its energy, which can in principle be obtained by solv- 
ing Eq. (O, determines the island charge difference at 
zero temperature, n(N, 0). 

For —Ej < A — 2Ec < the gap between the ground 
state and the bottom of the continuum is of order Ej. 
The gate voltage dependence of the ground state energy 
and the island charge remain qualitatively the same as 
those for A — 2Ec > 0. In particular, the charge plateaux 
with n = ±2 are joined by a smooth curve that lacks an 
intermediate plateau at n = 0. 

If A decreases even further, A — 2Ec < —Ej, the 
island charge difference develops a plateau, n(N,0) « 
for I AT I < 1 — A/2Ec- In this regime the ground state 
corresponds to a bound state of two quasiparticles, one in 
each island, with the gap between the ground state and 
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the bottom of the continuum being exponentially small, 
£ 11+ 2A-i5«Aexp (- Ig^fVl?-^)^ Nearthe 
charge degeneracy points N± = ± (l — je^) between 

the states with n = and n = ±2 Eq. 10 simplifies. As 
a result the island charge difference can be found from 
the implicit relation, 
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Here we assumed TV to be positive and denoted by 
N* = Z§f (N - 1 + ^j^J a rescaled gate voltage dis- 
tance from the charge degeneracy point. For N < 
the island charge difference can be obtained by noting 
that n is an odd function of N. Away from the charge 
degeneracy point, i.e. for \N*\ 3> 1, we obtain the ap- 
proximate expressions for the island charge difference: 

n = $f cxp f-li^j for N* < 0, and n = 2 - ^ for 
N* > 0. 

In summary, we have studied the low temperature ther- 
modynamic properties of a double-island qubit in the 
limit of the vanishing single particle mean level spacing. 
We have shown that the ground state of the system is al- 
ways separated from the continuum of excited states by 
a finite gap, and determined the gate voltage dependence 
of the island charge. 

For an odd number of electrons in the device the 
ground state corresponds to a bound state of the in- 
trinsic quasiparticle. The binding energy of this state 
is given by Eq. (@J. It attains its biggest value given 
by the Josephson energy Ej at the charge degeneracy 
point. The gate voltage dependence of the island charge 
is given by Eqs. © and ©■ The presence of the bound 
state leads to a non-monotonic temperature dependence 
of the width of the Coulomb blockade step. The minimal 



width is achieved at the ionization temperature of the 
bound quasiparticle state ~ Ej/{± In S7r f 2 Ej ). 

For an even number of electrons the gate voltage de- 
pendence of the ground state energy and the island 
charge difference are determined by Eqs. JUJ and l|12(l 
respectively. For A < 2E C and |JV| < 1 - A/2E C the 
two intrinsic quasiparticles are bound to the tunneling 
contact. The binding energy is exponentially small, see 
discussion above Eq. Ijl2|l . For A > 2Ec the Coulomb 
stimulation of the Josephson energy in the double-island 
system can significantly exceed that for a single grain 
coupled to a bulk electrode, see Eqs. IjlUII and l|llfl . This 
is in agreement with the experimental findings of Ref. |j] . 

In deriving our results we neglected virtual transi- 
tions to higher energy states, such as the non-resonant 
charge states and states with more than one quasipar- 
ticle in the system. Such transitions lead to corrections 
that remain small at the charge degeneracy point and 
amount to renormalization of the charging energy of the 
system [ijj, [l9j . Thus Eq should be understood as the 
renormalized charging energy. 

Our results can be directly tested by sensitive charge 
measurements techniques that became available re- 
cently 0, |2(|. The bound quasiparticle state is also 
present in a single superconducting grain coupled to a su- 
perconducting lead when the electrostatic energy favors 
an odd number of particles in the grain However in 
contrast to the latter case [lOj the gate voltage depen- 
dence of the island charge in the double-island qubit has 
no discontinuity at the charge degeneracy point. This 
difference is due to the absence of thermal quasiparticles 
in the double-island system at low temperatures. 
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